On the spectrum of facet crystallization waves at the smooth ^He crystal surface 
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The wavelike processes of crystallization and melting or crystallization waves are well known 
to exist at the *He crystal surface in the rough state. Much less is known about crystallization 
waves for the He crystal surface in the smooth well-faceted state below the roughening transition 
temperature. To meet the lack, we analyze here the spectrum of facet crystallization waves and its 
dependence upon the wavelength, perturbation amplitude, and the number of possible facet steps 
distributed somehow over the wavelength. All the distinctive features of facet crystallization waves 
from conventional waves at the rough surface result from a nonanalytic cusplike behavior in the 
angle dependence for the surface tension of smooth crystal facets. 

PACS numbers: 67.80.-s, 68.08.-p 
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I. INTRODUCTION 

Since the prediction by Andreev and Parshin 1] in 
1978, crystallization vifaves at the superfluid-solid ''He in- 
terface have become a well-knovifn phenomenon. At low 
temperatures the ^He crystal in contact with its liquid 
phase can support weakly damped oscillations of the in- 
terface due to processes of crystallization and melting. 
From the dynamical point of view such weakly damped 
crystallization waves at the rough crystal surface are 
an immediate counterpart of the familiar gravitational- 
capillary waves at the interface between two normal liq- 
uids (see review [2[). 

To date, the crystallization wave phenomena in ^He 
have extensively been studied for the rough state of crys- 
tal surfaces, but not much study has been made for 
the well-faceted and atomically smooth crystal surfaces 
which may have an infinitely large surface stiffness. The 
most distinctive feature of the smooth faceted surface 
from the rough one is the existence of nonanalytic cus- 
plike behavior in the angle dependence for the surface 
tension, e.g., P-i4l|. The presence of singularity leads to 
qualitative distinctions in a number of the phenomena at 
the smooth faceted crystal surface, e.g., amplitude de- 
pendence velocity of traveling waves 'S' , quantum finger- 
ing of the inverted liquid-crystal interface in the field of 
gravity [Q] , Ray leigh- Taylor instability with generation of 
crystallization waves [7[ . 

In the present work we develop a theory on the spec- 
trum of facet crystallization waves at the smooth faceted 
surface of a '*He crystal in contact with its liquid phase. 
We consider a few possible types of facet crystallization 
waves and determine the dispersion relation between the 
frequency and the wave vector, perturbation wave am- 
plitude and the number of the crystal facet steps per 
wavelength. For simplicity, we discuss the basal plane of 
hexagonal "'He crystal as an example of the crystal facet 
and disregard any anisotropy in the basal plane. 



II. LAGRANGIAN 

Let us assume that the crystal surface is parallel to the 
x-y plane, with vertical position at z = 0. In order to 
derive the oscillation spectrum of a facet surface, we pro- 
ceed as follows. Let us call ( — ({t, r) the displacement 
of the surface from its horizontal position z = 0. In the 
lack of energy dissipation the surface oscillations can be 
described by the action 



S^ J dtL[at,r),C{t,r)] 
with the following Lagrangian 



(1) 



L[C{t, r), C{t, r)] = 



Pet 



d^rd 



,2^, gt, r)C(t, rQ 
27r|r — r'\ 



d^r (^a(n) yiTTvcF + ^ Ap5C') ■ 
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Here r = (x, y) is a two-dimensional radius- vector. The 
first term in the Lagrangian represents the kinetic en- 
ergy of the interface having an effective density pd- We 
assume that both the liquid and the solid phases are in- 
compressible. Because of low temperature consideration 
we will also neglect the normal component density in the 
superfluid phase or, equivalently, difference between the 
superfluid density ps and the density of the liquid phase 
p. Then the effective interface density pof is given by 

Pet = {p' - Pf/P 

and depends on the difference Ap = p' — p between the 
solid density p' and the liquid density p. For our pur- 
poses, the exact magnitude of the effective density is 
inessential. 

The second and third terms in the Lagrangian are the 
surface energy and potential energy of the interface in 
the field of gravity with acceleration g. 

Unlike the liquid-liquid interface, the surface tension 
a{n) for the crystal facet depends essentially on the di- 
rection of the normal n to the interface. In our simple 
case this is a function of angle i? alone between the normal 



and the crystallographic [0001] or c-axis of the crystal hep 
structure with the geometric relation | tan 9 \ — \ VC |- 

For the crystal facet tilted by small angle d from the 
basal plane, the expansion of surface tension a{9) in se- 
ries in 6 starts as, e.g., Refs.0, |j]. 



a{6) = ao + ai tan | i? | + . 



tan( 



VC|. (2) 



We do not write the next terms of expansion, e.g., cubic 
one due to step-step interaction, since we assume to study 
only small bending of the crystal surface. The angular 



behavior has a nonanalytic cusplike behavior at 6* = due 
to «! = ai (T) representing a ratio of the linear facet step 
energy /3 to the crystallographic interplane spacing. Be- 
low the roughening transition temperature for the basal 
plane Tr '^1.2 K the facet step energy (3 — f3{T) is posi- 
tive and vanishes for temperatures T > Tr. 

To consider a traveling wave, e.g., propagating to from 
the left to the right, we represent the interface perturba- 
tion as C(i, r) = ({r — Vt). Here V is the phase velocity 
of the wave. Then the action ([TJ can be written as 



J 



-'-^f/i£^ 



>- ^ '-^^^) (^ ■ ^^^4^) - /^'.{"(")v/TT(W+ iA«c' 



r 



Since the integration is performed over r and r' within gunrent in ( by Vt. Next, by integrating twice by parts, 
the infinite limits and the kernel in the kinetic term de- we arrive at 
pends on difference | r — r' | alone, we can shift the ar- 



J 



S ■ 



dti^jd\d\'Pfar)ar'){v--^ 



dr 



27r|r 



-jd\ (a(n)^l + (VC)2 + \^pgC 



r 



In what follows, we will study sufficiently small bend- 
ing of the crystal surface with the sufficiently small dis- 
placements C and small angles 6. Involving inequality 
I VC 1^ 1 and I tan 6* | = | V^ |, we take only first terms 
in the expansion of the surface energy 



Next, we will choose the a;- axis as a direction of the wave 
propagation V = [V^ 0, 0) and replace x' with x in the 
spatial derivatives. Finally, we arrive at 



a(n)Vl + (VC)2« ao + ai I VC | +ao(VC)V2. 



S[Q~S[Q - 0] 



dt 



PctV' 



.jsrsr'am^')^. {.^hv)*IH-^ I '■■■ I +T<'«'+^4)}- 



(3) 



Variation 5S/5C,{x) yields the equation for interface os- 
cillations 



G(r - r')ax') d\' + a,^ («gn(£)) = « ' ^ 



Here, for convenience, we have introduced the Green 



function according to 



G{r-r') = -p,,V^ — 



1 



dx'^ V 27r|r — r' 
dx 



aQT^5{r - r') - Apg . (5) 



The solution of Eq. ^ ior V = has been studied in 
Ref.Q. 

Before solving Eq. Q, we make the following remarks. 
First, in the regions with ('{x) ^ the equation (|3]) 



reduces to a linear equation with the difference kernel 
G{r-r')(:{x')d^r' =0. 



The solution can be found as a sum of independent 
Fourier harmonics C{x) = X^n Cg 6xp(iga;). For a single 
harmonic 

C{x) = Cq exp{±iqx) or C{x) ^ Cq cos q{x - xq) , (6) 

one should have G{q)C,q — 0. For an existence of non- 
trivial solutions, it is necessary to put G{q) — 0. Thus, 
vector q must satisfy the condition 

G[q) = pcf ■^^'- - aog' - AP5 = Pciqiy - V^{q))= 0. 
q 

(7) 
Here we have introduced notation Vo(g) for the phase ve- 
locity of crystallization waves at the rough crystal surface 
with the spectrum ujQ{q): 



Vo{q) = UJo{q)/q , uJo{q) = \/{aoq^ + ^P9q)/p> 



cf 



Second, in the regions with C{x) = the solution is 
trivial, i.e., 

C(x) — const. 

And the last, since during the melting-crystallization pro- 
cess the total mass of the solid and liquid phases remains 
invariable, the solution of Eq. ((4]) must satisfy the follow- 
ing condition 



C{x) dx = 0. 



As a result, the general solution for profile C,{x) should 
represent a train of flat segments and half-sinusoids. The 
width of a half-sinusoid, which we denote 1/2, is governed 
by the magnitude of vector q according to G(g) = or 
Vo{q) = V and is equal to 1/2 — n/q. The half-sinusoid, 
which connects two neighboring flat segments, can be re- 
garded as a macroscopic facet step in contrast to elemen- 
tary steps of an atomic scale. To illustrate, we give two 
examples of such facet crystallization waves with alterna- 
tion of one or two various flat segments in Figs. [1] and [5] 
It is obvious that the wavelength A cannot be less than 
a sum of two half-sinusoids I = 2-k /q. In other words, the 
wave number fc = 27r/A does not exceed vector g, i.e., 
k ^ q. Later, it will be seen that always uj{k) ^ ujQ{k). 

Besides the wavelength, the wave spectrum at the crys- 
tal facet will also be characterized by the number of fiat 
segments and their amplitudes. The width of macro- 
scopic crystal steps 1/2 — ix /q depends on a ratio of wave 
frequency a; to wave number fc and can be determined 
using relation 



k 



'^o(g) 



-A/ 2 



A/2 



FIG. 1: The wave is formed by alternating one flat segment 
of length (A - l)/2 and kink of width 1/2. 
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FIG. 2: The wave is composed of alternating two flat segments 
of different lengths and four kinks of width 1/2 each. 



Provided I <C A, the width of half-sinusoids can be ne- 
glected in first approximations and the macroscopic facet 
steps can be treated as a kink of zero width. Usually, this 
is a range of sufficiently high frequencies uj{k) S> wo(A:) 
and small wave perturbation amplitudes | CI ^ ^ <C A. 



III. SIMPLEST TYPE OF WAVES 

Let us start from the simplest type of the crystalliza- 
tion wave which is formed by alternating a single fiat 
segment of height C and a half-sinusoid of width 1/2 so 
that / C{x) dx = (Fig. [TJ. In this case the wave pertur- 
bation is specified by 



x \< 



X-l 



Cix) = C <^ cos [2/L(| :, I _ A_i)] ^ A^ <| a; |< A±i . (9) 



1, 



Y<l^l<l 



(8) 



The total length of the fiat segments is equal to X — I. 
For correctness of approximation used, we should assume 
that I C |< ^271- =^ A/27r. 

To determine the spectrum for such type of the facet 
wave, we calculate the action as a function of the wave 
amplitude C and then minimize the action. The calcu- 
lation is readily performed in the Fourier representation 
using 

C(a;) = 5]Cne''="% 

n 

kn^2Trn/X^kn (n = 0, ±1, ±2 . . .). (10) 



Hence we have for the variation of the action taken from 
the flat crystal surface and ascribed per wavelength 

_. oo 
5[C]-5[C = 0] = - J2 G.(fcn)CnC-n-4ailC|- (H) 



with the obvious properties C-n = (n and Cn=2m = 0. For 
A = /, harmonics /i = /_i = 1/2 alone remain nonzero. 
Inserting Eq. (fT2|) into (fTTj) and minimizing the action 
dS/d I C 1= 0, we arrive at the equation which determines 
the spectrum of facet crystallization wave 



The second term is a contribution due to formation of the 
facet segment with the step amplitude 2(^. The number of 
such steps over wavelength equals 2 in our specific case. 
The calculation of Fourier components is straightforward 



A/2 



ax)e' 



■ dx = (Xfn , 



A/2 

sin(7m/2) / cos{Tml/2X) nn . 

^'' = (7rn/2) \l-nny\^ " ''°' T ' ^ ' 






(-'Lj{kn)Jn ^ 



2ai 



Then, using p^ and I — 2Tr/q, we have 



(13) 



^C-^^,- "o'^l^"' + ')' - ^«)/5."« - 5^ 



J2m+1 



aik 



1 



cos2[|(2m + l)^] 
7r2 (2m +1)2 [i_(2m+l)2fc2/^2- 

I 



2 ' 



m = 0, 1, 2, 



(14) 



The magnitude of vector q is determined from the con- 
dition ujo{q)/q = uj/k. 

We first analyze the limiting case of infinitely narrow 
step / <C A or fc ^ q. In this approximation one can 
here neglect the contributions from the regular ao surface 
and gravitational Apg terms and estimate /2m+i at point 
fc = 0. Then, a sum in (ITil) reduces to 



2 A °° 

uj 4 ^-^ 



1 



^°^ k TT-^ ^ (2m + 1)3 



703) 
27r2 



where C(3) « 1.20 is the Riemann zeta-function. Finally, 
we get the spectrum 



aifc2 



47r3 



0L\ 



7C(3)pcf|CI 7C(3)pefA2|C| 



(15) 



The phase velocity V depends on the perturbation am- 
plitude alone 



0L\ 



V=-=(— 

k V7C(3)pcf|C 



1/2 



A special feature of the spectrum is its growing stifi'ness 
as I C 1^- when the frequency of surface oscillations 
becomes infinitely large. Regardless of magnitudes oq 
and Apg the behavior w '^| C |~i/2 is universal in the 
I C K limit. 

To understand the typical magnitudes and applicabil- 
ity of our approximations, we start from the case when 



wavelength A exceeds slightly the double width of crys- 
tal facet step, i.e., (A - l)/X < 1. Then, in Eq. ([14]) the 
main contribution results from the first term with m — 0. 
Approximately, one has 



1 ai/c2 



1 



uj{k) = ojo{k) H — - — — 

TT PcfWo(fc) I C I 



(16) 



The effect of crystal facet smoothness on the wave spec- 
trum due to nonzero ai is governed by a ratio 



ai/c2 



1 



2TTai/aQ A 
aok^ + Apg Yla " 1 + A2/A2 K" 



where Aq = 2Try^ao/Apg is a usual capillary length. To 
have a strong effect on the spectrum in the sense u}{k) :s> 
a;o(fc), we should satisfy the inequality 



Cl« 



27rai 



ao 1 + A2/A2 



(17) 



Below, in detail, we will analyze most interesting re- 
gion of sufficiently small perturbation amplitudes C ~^ 
and wavelengths smaller than capillary length A ^ Aq. 
Neglecting gravitational term in ([13]), we have 



E 



1 



^^ V (2"^ + 1)^ 



cos2[f(2m-|-l) 



[l-(2m-Hl)2fc2/g2]" 
TT a\k'^ 



8 Pef ^2 I C I 



-A/2 




FIG. 3; The sketch of crystaUization wave of the total per- 
turbation amphtude (" and wavelength A. 



Involving that k/q — agk^ / (pcf ui'^) -C 1 and estimating 
the above sum within logarithmic accuracy as 

7C{3)/8~{7Ty8-l){k/qf\n{q/k), 

we find the spectrum with the correction due to finiteness 
of ao 



TT aik^ 
703) PcfIC 



21C(3) V 71- ao 

nai 1 



2/3 



X In 



7C(3)ao k\C 



The width of the kink between two flat segments can 
readily be estimated from q = pciUJ^{k)/{aQk'^) as 



1/2 « 7C(3) 






c 



(18) 



Thus, the approximation of zero-width facet step / ^ A 
can be justified for the small amplitude perturbations if 



Cl« 



A 



ai 



14C(3) ao 



On the whole, the spectrum of facet crystallization 
waves can qualitatively be described by introducing effec- 
tive surface tension or stiffness dependent on both wave 
vector and perturbation amplitude 



Otei 



ao 



7«'>lfr 



(19) 



IV. WAVES WITH A FEW CRYSTAL FACET 
STEPS 

Here we consider another type of facet crystallization 
waves with an arbitrary number of crystal facet steps per 
wavelength A. The wave amplitude C, is assumed to be 
sufficiently small in order to neglect the regular surface 
ao and gravitational /S.pg terms. This limit corresponds 
to high frequencies u){k) ^ ujo{k) and zero width of the 
kink between two flat segments. Thus, the profile of the 
perturbed crystal surface represents a broken line con- 
sisting of vertical steps and horizontal segments. For the 
definiteness, we consider the surface profile from AN ver- 
tical steps linking the same Af = AN flat facet segments 
with length Ax = X/M (Fig. [3]). The vertical steps, hav- 



ing the same height of C/N, are located at the points 

2fc — 1 

Xk = ±^— Ax, k^l,2,... ,N. 

The points Xk+i and Xk are connected with a horizontal 
segment of length Ax with the vertical amplitude 

N— \ k — N \ 
Ck = ±C '_^ ' , k^0,l,2,...,2N. 



In accordance with ()10|) we find the Fourier components 
for C(x) 

A/2 

Cn= /"dxC(a;)e-2™^/^ 

-A/2 

2Ar-l "''=+1 

E N-\k-Nl / , , 

C « / sin(27rnx/A) dx 



k=l 



N 



= "*— wZ^(''°'^^"^°''2/fe+iv), yfc = — -(2fc-l). 



fe=i 



Summation, in essence, is reduced to summing geometric 
series. As a result, we obtain 

2Trn N sin[nn/4:N) 

with the obvious properties C-n — Cn and Cn=2m — 0. 

Unlike derivation of Eq. (fTTj) . we must take here into 
account that each crystal step contributes ai | C | /N 
into the action and the number of steps equals AN at the 
wavelength. Since ANai \ C \/N — Aa\ \ C |, we have 
the same contribution to the action from the steps and 
thus the same form ([T5| of the equation to determine the 
wave spectrum 



E 

m— 



/OofW^ I /2m+l ? Oiik 



k (2TO+1) 27r I CI 



(21) 



but with another Fourier component /2m+i- Applying 
/2m+i from (pn|) . we have 



oo _. 

V L_ 

^^ (2m + 1) 



1 



TT aifc^ 



^^g (2m + 1)3 N^ sin2[7r(2m + 1)/AN] 2 pef |C 1' 

Finally, the spectrum of facet crystallization waves is 
determined by 



2 27r aifc 



7C(3) Pcf I C I ' 



(22) 



where sjy is given by 



4 J2 (2m +1) 

m=0 



'AT 



N^ Y, (2m+l)-3sin"2[7r(2m+l)/4A^] 

m=0 



TABLE I: The ratio of frequencies sn ~ un /'^i for facet 
crystallization waves with the different number A^ of crystal 
facet steps per wavelength 



A^ 


1 


2 


3 


4 


5 


6 


7 


oo 


Sn 


1.000 


1.103 


1.122 


1.128 


1.131 


1.133 


1.134 


1.136 



and 

"""^ 62C(5) ■ 

The factor sj^ varies insignificantly within the range from 
1 for A^ = 1 to 1.136 for iV = oo (TableU]). 

We have analyzed above the wave spectrum for the 
regular arrangement and identical height of crystal steps 
over wavelength. This is not, of course, solely possible 
structure with AN crystal steps. The location of crystal 
steps and their heights can have an arbitrary and disor- 
dered structure. However, the dimensional estimate ((22|) 
for the spectrum holds for. As concerns the factor sn, it 
varies slightly as a function of the perturbation profile. 

Provided the regular part of surface tension ao differs 
from zero, the maximum number of crystal steps 4iVm 
is limited. Using magnitude I for the width of the kink 
between two flat segments 



I 






CI, 



we estimate the number ANm of possible crystal steps for 
the given frequency w according to 



ANr. 



A 

T 






(23) 



Obviously, the smaller the perturbation amplitude | C l> 
the larger the number of possible crystal steps. 

On the whole, the frequency of crystallization waves at 
a smooth crystal facet proves to be dependent not only on 
the wavelength, but also on the wave amplitude and the 
number of crystal facet steps which can be placed within 
the wavelength. From the experimental point of view this 
means that the excitation of crystallization waves with a 
fixed frequency will result in exciting some train of waves 
with different wavelengths, amplitude, and the number of 
crystal facet steps. In this connection the shape of the 
perturbed crystal facet will resemble rather an irregular 
and ill-defined profile with some elements of irregular- 
like character. A weak dependence of the wave frequency 
on the number of steps facilitates such phenomenon. In 
some sense one might say about transition to a rough 
state of the surface and destruction of the crystal faceting 
0- 



CONCLUSION 



cated picture than that at the rough crystal surfaces. 
The plane crystallization wave represents an alternation 
of flat crystal facets linked via macroscopic crystal steps 
with the width dependent on the wave velocity. Most 
striking phenomena should appear in the limit of suf- 
ficiently small perturbation amplitudes. The frequency 
spectrum of facet waves depends significantly on the per- 
turbation amplitude. The dependence on the structure 
of wave perturbation and the number of crystal facet 
steps is not so drastic. Excitation of waves at a given 
frequency should produce a train consisting of waves dif- 
fering in wavelength, structure and the number of facet 
steps and interacting nonlinearly with each other. This 
all is in contrast to harmonic waves which exist at the 
rough crystal surfaces and have the lower frequencies at 
the same wavelengths. Evidently, the distinction results 
from a singularity in the angular behavior of the surface 
tension for the smooth crystal facets. 

Let us estimate typical frequencies for the short wave- 
length range A ^ Aq '^ 6 mm. Taking pof ^ 2 mg/cm"^, 
ai ~ 0.014 erg/cm^ and ao ~ 0.16-0.18 erg/cm^, e.g., 

for the (0001) '*IIe facet, we have the frequency uj ^ 
10 kHz and velocity 1/ ~ 2 m/s for the wavelength A ~ 

1 mm and perturbation amplitude | C |~ 1 /.tm. In this 
case one may expect the maximum number of possible 
steps over wavelength to ^ 100. If the perturbation am- 
plitude for the same A approaches | C l~ 0-1 nrm, only 
one or two steps become possible. The wave frequency 
reduces to about w ~ 2 kHz which insignificantly ex- 
ceeds the magnitude at the rough surface. Note that, for 
the perturbation amplitude of a crystal lattice spacing 
in height, the propagation velocity V reaches the magni- 
tudes of '-^lOO m/s comparable with the sound velocity. 

In principle, one can find a few casual mentions about 
phenomena similar to the formation of crystallization 
waves at the crystal ^He facets under heavy shake of an 
experimental cell [8| or in the process of anomalously 
fast growth of a '^He crystal under high overpressures 
[3, [ifl]- More convincing observation in favor of an 
existence of progressive facet waves has recently been 
found |11,] at the crystal (001) facet in ^He. Apparently, 
one of complicating factors in exciting and studying 
facet crystallization waves may be associated with the 
threshold character for most of phenomena occurring at 
the smooth crystal facets. In particular, it may require 
a sufficiently large size of the facet and sufficiently high 
amplitudes of driving perturbation. In this connection 
it may be helpful to employ the conditions close to an 
onset of some instability, e.g., electrocapillary one in an 
electric field across the interface^, Ray leigh- Taylor Q 
or Faraday instabilities [12 |. 



The crystallization waves at the smooth crystal facets 
are expected to demonstrate a more varied and compli- 
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